Phonon spectroscopy through the electronic density of states in graphene 
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We study how phonon structure manifests itself in the electronic density of states of graphene. 
A procedure for extracting the value of the electron-phonon renormalization A is developed. In 
addition, we identify direct phonon structures. With increasing doping, these structures, along with 
A, grow in amplitude and no longer display particle-hole symmetry. 
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It is a remarkable result of many body physics that in 
electronic systems in which the density of states (DOS) is 
constant on the scale of a phonon energv electron-phonon 
renormalizations entirely drop out[l|,l2| and no phonon 
signatures are expected or seen. This has now changed 
with the advent of graphene. Graphene was isolated only 
m 20040, H but has since been extensively studied and 
found to exhibit many unusual properties. These include 
a novel quantum Hall effect, a minimum conductivity, a 
Berry phase of tt and other effects related to the chirality 
of its charge carriers. 0, 0, 0] 

In this letter, we show how measurements of the 
DOS in graphene offer a new opportunity to obtain de- 
tailed information on electron-phonon coupling in sharp 
contrast to the case of ordinary metals. This is im- 
portant as it illustrates a strong violation of a well- 
established result of many body physics and provides 
unusual doping-dependent predictions for the manifes- 
tation of the electron-phonon interaction in experiments 
on graphene. This result arises for two reasons. First, 
the charge carriers exhibit relativistic dispersions with 
quasiparticle energy (ek) linear, rather than quadratic, in 
momentum (k), ek — ±?iwo|k|. Here, plays the role of 
the velocity of light and the ± gives the upper and lower 
Dirac cones, respectively. At neutrality, the lower cone is 
fully occupied and the upper one, empty. This dispersion 
gives rise to an energy dependence of the DOS which is 
linear to high energies. Second, both theory and experi- 
ment indicate that the major coupling is to high energy 
phonons of order 200 meV.0, S Hoj Thus, an electron 
scattering from an initial state to a final state through the 
assistance of a phonon will sample changes in the DOS 
on the scale of the phonon energy, which is significant 
in graphene. Another special feature of graphene is that 
the number of charge carriers can be changed through 
charging in a field effect device where the chemical po- 
tential (/io)i measured with respect to the Dirac point 
(neutrality point), is proportional to the square root of 
the gate voltage. These characteristics offer a rich new 
spectroscopy for the study of phonon effects including 
their variation with /ip- 

In graphene, the bare band DOS varies linearly with 



energy and hence the renormalized density of states N{lo) 
is given by [HI, [13, El, E | 



N(uj) 



Wc 



dt'- 



-ImE(a;) 



_wc ^ Rel](w) -f ^ - e]2 + [ImE(t^)]2 ' 

(1) 

with iVo = In Eq. Ilj, Wc is an upper cutoff 

on the Dirac cones given by yiv^t, with t the nearest 
neighbor hopping parameter, and S(a;) is the electronic 
self-energy given by[ll| 
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where a^F{v) is the electron-phonon spectral density and 
n{v) and /(w') are, respectively, the Bose-Einstein and 
Fermi-Dirac distribution functions at temperature T. For 
the bare band case (S([x)) 0), the chemical potential /i 
reduces to its noninteracting value /io and sets the doping 
level. Also, the Lorentzian form in Eq. ([T|) reduces to 
5(lj + fiQ — e) and the DOS becomes \uj + fio\- 

Park et al.[8] performed a full first-principles study of 
the electron-phonon interaction in graphene and found 
that the result could be approximated by an Einstein 
mode at 200 meV. For coupling to an Einstein mode, 
the electronic self-energy in a system with a linear DOS 
can be evaluated in the usual manner to give, at zero 
temperature, an analytic form with 15 1 
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where A is the area under the Einstein mode and lue is 
the Einstein frequency. For simplicity, we have assumed 
in writing Eq. ([3]) that Wc is larger than any other energy 
of interest, but in all numerical results presented here, 
this approximation was not made. The corresponding 
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imaginary part is — ImS(w) 

io < Wc - fJ.0 
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uj > —Wc — fiQ — ^E- In terms of this self-energy, the 
renormaUzed density of states is given by 
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where F = —ImT,{uj) and uj = uj — ReE(a;) -t- fi. For 
finite /i, the problem no longer has particle-hole sym- 
metry and Rel](a; = 0) is not zero and provides a shift 
in chemical potential from bare to dressed value with 
M = Mo + Rei;(a; = 0).[l3, For a clean system F will 
vanish for ~uje < uj < loe, and Eq. ([4]) reduces to 



N{uj) 



— ujsgiuj, for — lue < lu < lue- 



(5) 



In this special range, the DOS is very closely related to 
the ReE(w). Returning to Eq. ([T]), it is important to re- 
alize that, for infinite bands with constant DOS, the |e| 
factor would not appear and the integral over e would 
give a constant independent of to so that phonon renor- 
malizations simply drop out. Graphene is very different. 
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FIG. 1: (Color online) —dYieTiiu)) / duj vs ui for a truncated 
Lorentzian electron-phonon spectral density peaked around 
uje ~ 200 meV. Shown are curves for = 0, 150, 500, and 
700 meV. 

While Eq. ^ has been written for a single Einstein 
oscillator, it nevertheless provides us with valuable in- 
sight into the relationship between phonon structure and 
the real part of the self-energy. Rel](ijj) has singulari- 
ties of the form Injo; ± at uj ~ ±uje and a third 
weaker logarithmic singularity of the type {uj + uje + 
/io) In + u! + uje\ at ciJ = — (/io + loe)- The neutrality 
point is special, however. For — Q, only two singu- 
larities remain and they are both of the weaker kind 
{lo ± In |w ± uje\- In a real system there will of 



course be a distribution of phonons and the self-energy 
of Eq. ^ needs to be averaged over such a distribu- 
tion. This will reduce the prominence of the expected 
singularities in this quantity. In such a case it becomes 
useful to consider a first derivative — dReS(a')/du;. This 
is shown in Fig. [T] for four values of the chemical poten- 
tial, /io = 0, 150, 500 and 700 meV. The phonon distribu- 
tion used in these numerical calculations was a truncated 
Lorentzian centered around uje = 200 meV with width 
(5 = 15 me V. [13, [13] ■ As expected the top left frame ex- 
hibits only two phonon anomalies while the three other 
frames have three. Also in these three cases the anoma- 
lies at cij = zbcjB are much more pronounced than the 
ones at lo — ~[uje + ^lo), and also than those in the top 
left frame. In all four frames, the black dotted horizontal 
line was drawn through the local minimum at w = and 
identifies the value of the electron-phonon mass renor- 
malization parameter A as we will now describe. For u) 
small near the Fermi energy {uj = 0), ReS(a;) in Eq. ([3]) 
can be shown to vary as ReE(w) ~ —\uj -\- ReE(w = 0) 
and the dressed quasiparticle energy i?k is given by the 
equation i^k — ReE(i?k) + /i = ±?ifo|k| = i?k(l + A) -|-/io. 
OvE^,= [±hvo\k\-^io\l{l + X) = ±nvo{k-kF)l{'i + X), 
which means that A simply renormalizes the bare Fermi 
velocity from tiQ to Vq = vo/{l + A). As Fig. [1] shows, 
A grows with increasing /zq as Eq. ([3]) implies. The red 
arrow indicates the Dirac point defined by |k| = 0. 
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FIG. 2: (Color online) (a) iV(a;)/iVo (in eV) vs uj for the 
bare chemical potential fio = 150 meV. The solid blue curve 
gives the phonon renormalized case and the black dotted, the 
bare band case. The arrows show the bare (long black) and 
renormalized (short blue) value of jj,. (b) Same as for (a) but 
with /io — 500 meV. 

The relationship between boson structure in the self- 
energy and its manifestation in the DOS is given by 
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Eq. (g]). Results for N{ll:) are shown in Fig. [2l The 
frame (a) is for fiQ — 150 meV, which is smaher than loe 
and (b) is for fiQ — 500 meV > uje- The shaded yel- 
\ow region is the occupied part of the bare band which is 
shown as the black dotted curve. Phonon renormaliza- 
tions change the shape of the DOS and hence the value of 
the chemical potential must be altered to keep the correct 
number of particles. The long black and short blue ar- 
rows point to the value of the bare and dressed chemical 
potential with fi — fio — Rel](w = 0). Phonon anomalies 
are clearly seen in the dressed curves. To emphasize this 
structure an Einstein spectrum was used with cue — 200 
meV so that the phonon structures fall at a; = iiOE, one 
on either side of the Fermi energy. The expected singu- 
larity at cj = —{coe + fJ-o) is by comparison very weak and 
appears as a slight change in slope in Fig. (5] Two addi- 
tional features of these curves are to be noted. At the 
Fermi energy {uj = 0), the dressed and bare DOS have 
exactly the same value. In the region of the Fermi energy 
Eq. ^ applies and N{uj)/No = \lu{1 -|- A) -I- ^o|, which 
differs from its bare value only by the additional factor of 
(1 -I- A). At w = 0, this difference disappears and dressed 
and bare DOS are the same. Phonons do not change the 
value of the DOS at the Fermi level. The slope out of 
u! = 0, however, is changed by a factor of (1 -I- A) as can 
be seen in both frames of Fig. [2] and we also note that 
this linear behavior persists over a considerable energy 
range set by the value of the Einstein oscillator. Recog- 
nizing that the normalization for the DOS is Nq ~ I/wqj 
one might naively think that the (1 -f A) renormalization 
can be included in N{uj) simply by changing vq to Vq in 
Nq, but we see here that this is not correct. Only one 
(1-1- A) factor enters and not its square. The basic reason 
underlying this fact is that the coherent part of the elec- 
tronic Green's function, which defines the quasiparticles 
in the interacting system, contains only 1/(1 + A) of the 
spectral weight. The remainder A/ (1 + A) is found in the 
incoherent piece describing phonon- assisted processes. 

Phonon structure in N{uj) can be brought out through 
differentiation. Results for dN(uj)/duj vs uj are given in 
Fig. [3]as the solid blue curves, where N{lu) is normalized 
by Nq. Frame (a) is for /io — 150 meV and (b) is for 500 
meV. The vertical drop where dN (lo) / dio goes from pos- 
itive to negative is at the Dirac point of the interacting 
system. Comparison with the bare band case, the dotted 
black line, shows a small shift of the position of the Dirac 
point between bare and dressed case. The bare case pro- 
vides a useful reference line about which the effects of 
the electron-phonon interaction are easily seen. Besides 
the phonon structures at a; = ±lje, — (wb -I- /io)j we note 
that the height of the curve above one at uj — gives 
the value of A directly which increases significantly with 
increasing value of chemical potential [as shown by the 
red solid curve in (c)]. The red dashed fine is included 
for comparison and gives —dKeT,(Lo)/dLO. There are some 
differences between these two sets of results but we can 
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FIG. 3: (Color online) (a) dN{uo)/du vs uj (solid blue curve) 
for /io = 150 meV. The black dotted curve sets a baseline and 
is the bare band case. The red dashed, which is for compari- 
son, is [1 — dReE(aj)/d(x;]sgn(aj ~\-ujd)- (b) Same as for (a) but 
with ^0 = 500 meV. (c) The absolute height of the phonon 
peaks at about uj = loe (solid blue dots) and — (/lo + uje) 
(open red circles) along with the A variation with /io, the lat- 
ter indicated as the red solid curve using the righthand axis 



conclude that all qualitative features seen in the DOS 
curves can be seen in the ReI](Lj). This is not to say that 
the imaginary part of S(ct;) plays no significant role. In 
frame (b), we see clearly that the jump at the Dirac point 
energy is no longer vertical but exhibits some smearing. 
This can be traced to the behavior about the Dirac point 
in the DOS shown in Fig. [DJb). The DOS no longer 
goes to zero at this point {co^i) but rather has a minimum 
about which it rises as a quadratic {uj — ujdf', seen in 
experiment ,l3|- We can show that for \{uj — ujd)Z\ <^ F, 
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No ^~ ^ 
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with Z 



[dRe^{co)/dco]^=^ 



and F 



— [ImE((jj)](j=j^^, which shows the lifting of the Dirac 
point and its conversion from linear to quadratic in 
{uj—ujd). This immediately leads to the smearing at Dirac 
point noted in the blue curve of Fig.[3Ub). In Fig.[3ljc), we 
plot the absolute value of the height of the phonon struc- 
tures as a function of /ip for the a'^F{v) spectrum used 
here, a truncated Lorentzian (see in the inset of Fig. HI 
long-dashed red curve). While the height of the phonon 
peak at uo = —(uje + IJ-q) hardly changes with doping 
(/io) the other two peaks do, note the curve for uj = uje- 
We have not plotted the peak height for uj = —uje as 
for fiQ < CUE it is similar to the result for co = uje and 
for fiQ > uje it becomes ambiguous. These predictions 
provide verifiable tests that observed structures are in- 
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deed due to phonons. They also show how the increase 
in the DOS at the Fermi surface with increasing doping 
is reflected in larger coupling to the phonons. 
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FIG. 4: (Color online) (a) 1 - \dN{uj) /duj\ vs cu for = and 
uje ~ 155 meV (solid blue curve). The inset on the lower right 
compares the phonon region with the input electron-phonon 
spectral density (red-dashed curve), (b) 1 - ^N{uj)/No{lo) 
vs uj (solid blue curve) compared with the result using the 
procedure of Li et al. (long-dashed red curve) . 

In Fig.m^a), we show results for {l—\dN{uj)/du!\) in the 
specific case of /uq = which is close to the case recently 
observed in the scanning tunneling miscroscopy (STM) 
results of Li et al.Q To conform with those experiments, 
we have used ue = 155 meV, with lv = ±u>e shown as the 
vertical black dashed lines. The horizontal black dotted 
line through the local maximum of the solid blue curve 
at = identifies the value of A which we took to be 
0.2. A second important feature is the phonon structure 
which reflects the underlying a^F(i^) spectrum used. It 
is shown in the inset as the long dashed red curve where 
it is scaled down and compared with the absolute value 
of the blue curve about a; = 155 meV. While there is 
some agreement, the two curves have different profiles 
with the blue solid one much broader than the red long- 
dashed one. It is clear that such a plot is very useful 
in identifying phonon structure, i.e., not just the value 
of the mass enhancement factor A involved but also the 
position of the peaks in a'^F{i') and their strength. In ex- 
periments, it may be more desirable not to differentiate. 
In Fig. [IJb), we show as the solid blue curve a differ- 
ent quantity 1 — y/N{uj)/No{uj), where No{uj) is the bare 
band density of states. In this quantity, the value of the 
local maximum at w = gives 1 — vT+A — —0.095 
rather than the —A of frame (a). We also note that 
the phonon structures at oj = ±uje are not as sharp, 
however, some signature of a sharp peak in the a^F(v) 
used remains. The long-dashed red curve is for compar- 
ison and represents the quantity that was used by Li et 
al.[§| in their analysis of their STM data. They use a 
definition of an effective Fermi velocity dependent on E 
based on an integration of their conductance. They de- 
fine vf = dE/hdk with fc = ±| /j^ N{uj)duj\^/'^ , where 
is the energy locating the Dirac point. Like the blue solid 
curve, the maximum at w = provides 1 — y/1 + A and 
the phonon structures at w = ±u!e are clearly seen. The 



Li et al.0 estimate of A = 0.26 is close to the 0.3 value 
from ARPES.[i3| Our own estimate based on the STM 
data is somewhat higher, but carries considerable uncer- 
tainty because of the experimental noise. An important 
point to note between Fig. HJa) and (b) is to reiterate 
that the many body renormalizations correct N{lu) by a 
(1 -I- A) factor. Assuming \e\/vp l^l/Vf^ would over 
estimate the correction by an additional factor of (l-fA). 

In contrast to the standard expectation in wide band 
metals with nearly constant DOS on the phonon en- 
ergy scale, phonon structure does appear prominently 
in the DOS of graphene and this can be used as a new 
spectroscopy for determining electron-phonon coupling. 
We find that the mass enhancement parameter A can 
be extracted directly from the data around the Fermi 
energy. In addition, three prominent peaks are identi- 
fied ai Lu = ±uje and uj — ^{lue + fJ-o) associated with 
each Einstein mode. The size of these additional struc- 
tures increases with increasing doping as does the mass 
enhancement A. These increases reflect the increase in 
the underlying DOS. In this sense, graphene provides a 
new laboratory in which to study variations in electron- 
phonon coupling with changing carrier concentration. 
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